Surface charge relaxation and the pearling instability of charged surfactant tubes 
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The pearling instability of bilayer surfactant tubes was recently observed during the collapse of 
fluid monolayers of binary mixtures of DMPC— POPG and DPPC— POPG surfactants. We sug- 
gested it has the same physics as the well-known Raleigh instability under the action of the bilayer 
surface tension whose magnitude is dictated by the electrostatic interaction between charged sur- 
factants. In this paper, we calculate the relaxation of charge molecules during the deformation of 
the tubes into pearling structure. We find the functional dependence of the relaxation energy on 
the screening length explicitly. Relaxation effect lowers the cost of bending a tube into pearls 
making the cylindrical tube even more unstable. It is known that for weak screening case where the 
tube radius is smaller than the screening length of the solution, this relaxation effect is important. 
However, for the case of strong screening it is negligible. For the experiments mentioned, the situa- 
tion is marginal. In this case, we show this relaxation effect remains small. It gives less than 20% 
contribution to the total electrostatic energy. 

PACS numbers: 68.10.-m,61.30.-v,82.70.-y,87.22.Bt,02.40.-k,47.20.-k 



I. INTRODUCTION 

The formation of surfactant tubes and budding of 
spheroidal structures are of significant interest in biolog- 
ical processes. In particular, such structures constitute 
intermediates that are responsible for critical cellular pro- 
cesses like material trafficking from the Golgi complex 1 , 
and fusion and fission of membranes 2 . As seen during cell 
locomotion and the formation of Golgi structures, natu- 
ral surfactant tubes are prone to transform to a structure 
resembling a string of pearls. 

Pearling has been induced in tubular phospholipid 
membranes by adsorption of oifi or polymer—, on the one 
side of the membranes. These phenomenon were inter- 
preted in terms of the creation of membrane spontaneous 
curvature due to those external stimulus. 

We have recently observed pearling in tubular struc- 
tures formed during the collapse (2D-3D transition) 
of fluid monolayers of mixed phospholipids 5 . Col- 
lapse in binary monolayers of 7DPPC:3POPG and 
7DMPC:3POPG lead to the formation of cylindrical 
tubes 6 . These tubes can be 10s of microns in length, 
with diameters close to 1/mi (limit of resolution). A 
few of these are wide enough to resolve detailed features. 
As seen in Fig. such tubes show instability towards 
pearling without the introduction of any external gra- 
dients that may affect or induce the spontaneous curva- 
ture. Furthermore, the tubes, being microscopic and sub- 
merged in water are likely to be composed of surfactant 
bilayers, which arc in the liquid phase at the tempera- 
ture measured. This suggests that the tube surface does 
not have intrinsic spontaneous curvature itself. Thus the 
above mentioned mechanisms of pearling instability is 
questionable for the present case. 

In the same paper, we proposed that the pearling insta- 
bility observed is due to a simpler mechanism. Namely, 
the instability is caused by the surface tension energy of 
the surfactant bilayer. This is very similar to the well- 



FIG. 1: Two snapshots taken within a few seconds of each 
other, showing a surfactant tube undergoing a pearling insta- 
bility. The monolayer are 7DMPC3POPG binary mixture. 



known Raleigh instability of cylinder of fluid. Indeed, 
because one type of surfactants used in experiments are 
POPG surfactant, a charged molecule, the surface ten- 
sion of the bilayer should be, at least, of the same order 
of magnitude as the electrostatic energy per unit area 
of the tube. The latter is 7t<7 2 /kD where Co is the sur- 
face charge density of the bilayer, k is the inverse screen- 
ing radius of the solution and D = 80 is the dielectric 
constant of water. Using relevant experimental param- 
eters, this electrostatic energy is estimated to be about 
10 rriN/m. On the other hand the bending energy of a 
surfactant phospholipid bilayer is knownia to be about 
7 ~ 30/cbT (here, ks is the Boltzmann constant and T is 
the temperature of the solution). For a tube with radius 
Ro of 1 micron, this translates into an elastic (bending) 
energy of about 7/iZ§ — 10~ rnN/m per unit area. Thus 
the elastic energy is more or less negligible in comparison 
to surface tension energy. In other words, an instability 
similar to Raleigh instability of a fluid cylinder must be 
present for micron size tubes. For the case of a fluid 
cylinder, this instability leads to the breaking up of the 
cylinder into small droplets. However, for a surfactant 
tube, the breaking process is improbable because all sur- 
face tension energies involved are far below the rupture 
tensile stress (about lmN/m) of the lipid bilayer. The 
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pearled structure is obviously the most likely candidate 
for the final structure of this instability. 

We note that electrostatic surface tension energy is a 
negative quantity (repulsions between charged molecules 
prefer to expand, not to reduce, the surface area of the 
tube). Therefore, it cannot be the driving force of insta- 
bility. However, because the total surface tension of the 
bilayer must be positive, the presence of the negative sur- 
face tension forces the non-electrostatic surface tension to 
be positive and greater than the electrostatic one in mag- 
nitude. Quantitative calculation of the non-electrostatic 
surface tension from the balance between these energies 
requires the knowledge of the lateral compressibility of 
the bilayer and is presented in details in Ref. 0. It is 
shown that this non-electrostatic surface tension is of the 
same order of magnitude as the the electrostatic one, and 
thus, is much bigger than the bending energy. It is this 
non-electrostatic surface tension whose magnitude is dic- 
tated by the electrostatic counterpart which drives the 
tube towards pearling instability. 

This explanation of pearling instability as a Raleigh 
instability due to electrostatic effect is supported further 
by the experimental observation that adding monova- 
lent salt to the solution diminishes electrostatic energy 
and remove the instability of the cylinder. It is also ob- 
served that cylinders with smaller radius are more stable 
against pearling. This can be explained simply by the 
rapid increase of the elastic energy when the tube radius 
decreases. 

There is another well known electrostatics-induced 
pearling instability in literature, namely, the pearls-on- 
string structure of polyelectrolyte in poor solvent or of 
polyampholytes^. The physics behind these instabilities 
is different from that described above. Even though both 
instabilities are induced by electrostatic effects, in sur- 
face tension induced pearling, the characteristic size of 
pearls is determined kinetically. On the other hand, the 
pearl size and period of polyelectrolytes is determined 
thermodynamically by the balance between electrostatics 
energy and non-electrostatic energy (entropic or solvent- 
monomers interactions). This leads to the strong redis- 
tribution of the charge molecules (strong charge relax- 
ation) in the system to lower its overall free energy, which 
is obviously not needed in the case of Raleigh instabil- 
ity. One expects charge relaxation is the driving force of 
pearling instability when screening of the solution is weak 
such that the Debye screening length is much larger than 
the pearl size. On the other hand when the screening is 
strong, the Debye screening length is much smaller than 
the pearl size, electrostatics is a short range interaction 

and the instability is of the dynamical Raleigh type. 

n 

In the experimental system of Ref. .5, pure water is 
used. In this case, the screening radius of the system is 
comparable to the tube radius (about one micron). Thus, 
the situation is marginal, and it is not clear whether or 
not the relaxation of charge surfactants still plays a sig- 
nificant role in the pearling instability. In this paper, 
we would like to address this question by calculating ex- 



plicitly the gain in the electrostatic energy of the system 
when the charged molecules of the bilayer redistribute 
themselves during the pearling transition. We show that 
in this marginal case where the screening radius is equal 
to or smaller than the tube radius, the charge relaxations 
remains small. The energy gain due to this effect con- 
tributes at most 20% to the total electrostatic energy. 
For smaller screening radius, the ratio between these two 
energies decreases very fast (as fourth power in the ratio 
between the screening length and the pearl size). Thus 
in the experiments of Ref. 0, one can neglect the mod- 
ulation in surface charge density when considering the 
electrostatics of the system. 

This paper is organized as follows. In Sec. II, us- 
ing linear analysis, we briefly calculate changes in the 
surface tension and bending energies when a cylinder de- 
forms into a string-of-pearls. In Sec. Ill, we calculate the 
change in the electrostatic energy under this deformation 
and separate the contribution due to the relaxation of 
charged molecules. The latter is always negative. This 
gives an additional gain in the energy of deformation, 
making the tube even more unstable. In section IV, we 
discuss the relative importance of charge relaxation effect 
as well as various approximations involved. 

II. ELASTIC ENERGY 

Even though the elastic (surface tension and bending) 
energy changes when a cylinder undergoes pearling de- 
formation has been calculated^, before, we briefly repeat 
the calculation here in order to introduce the notations 
and to simplify their comparison with electrostatic en- 
ergy in later sections. 

Let us start with a model elastic free energy describing 
the cylindrical tubes. Denoting the bilayer tube length 
L, area S, and volume V, our starting free energy is the 
sum of the surface tension energy, the bending energy 
and an osmotic pressure energy: 

E = E s + E b + E =a J dS + 

+ J dS[2-/H 2 + *fK] +SpJdV (1) 

where 7 and 7 are the bending rigidity and the Gaussian 
bending rigidity of the bilayer, H and K are the mean 
and Gaussian curvature of the tube surface, dp is an os- 
motic pressure difference between the inner and outer 
region of the tube. In the above model, the osmotic pres- 
sure term is somewhat artificial. This term is needed to 
make the cylindrical shape the minimum of the energy for 
certain range of the parameters a, 7 and Sp (because the 
growth of the tubes is slow in experiments, we consider 
the tubes are in (quasi-) equilibrium and their shapes are 
determined by the minimum of the free energy). With- 
out this term, a spherical vesicle will always be the shape 
which minimizes the free energy, Eq. Q). This osmotic 
pressure term was also used by the authors of Ref. |9| to 
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study instability of cylindrical vesicles. A second choice 
for the model energy is to replace the osmotic pressure 
term in Eq. Q by a line tension term, which has been 
used by the authors of Ref. to study coiling instability 
in multilamellar tubes. Each of these models incorporate 
different physics in stabilizing the cylindrical tube. The 
choice of one model over the other is not important in 
this paper because we do not allow either the volume or 
the length of the cylinder to change in our subsequent 
analysis of the cylinder instability. In the free energy, 
Eq. 0}, the surface tension and the osmotic pressure are 
actually Lagrangian multipliers which enforce the restric- 
tions of area and volume conservation of the surfactant 
tube. 

Within linear analysis, to investigate the change in the 
energy of a tube undergoing a pearling instability, let 
us slightly deform the cylinder radially with a relative 
amplitude, e« 1, and a wave vector k (see Fig. As 
a result, the radius of the new tube varies along its axis 
according to: 



R(z) = R[l + ecos(kz) 



2n/k 



(2) 




FIG. 2: A cylindrical tube is deformed radially with wavevec- 
tor k. 



Due to the small permeability of water through the 
bilayer surface, in our model, we require that the tube 
volume does not change under deformation. This makes 
the average radius, R. of the new tube different from 
the original radius Rq. This new radius can be easily 
calculated. The volume per unit length of the new tube 
is: 



7 rir/k l-R(z) 

V = — / dz r> — J d2/1 1 - 2 

2tt J-Tr/k Jo 



2-nrdr = ttR 2 (1 + e 2 /2) (3) 



The condition of volume conservation then leads to the 
simple relation: 



R = R a /^l + £ 2 /2~ R Q (l - e 2 /4) 



(4) 



Let us calculate the change in the surface tension en- 
ergy of the tube. The element of the area of the new tube 
is, 



dS(z, 4>) = dzd4> R{z) v / l + R'{z) 2 



= dzd(f> R[l + e cos(fcz)] Jl+ [ekR sin(fcz)] 2 (5) 



where R'{z) = dR/dz. The area of the deformed tube 
per unit length is then, 



.4 
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2^ 



2vri? 



2vri? n 



irfk p2n 

/ dS{z,<j>) 

-n/k JO 

E[-(kR Q ) 2 e 2 /(l + e 2 /2)] 



(7r/2)Vl + e 2 /2 
1 - (kRo) 2 



for e < 1. (6) 



where E is the complete elliptic integral of the second 
kind. The change in surface tension energy per unit 
length is easily calculated to be: 



AE S = a(A-2irR ) 



TraRo 



[(kR ) 2 - l]e 2 



(7) 



Let us next calculate the change in the bending en- 
ergy. Because the Gaussian curvature energy 7 J dAK is 
a topological invariant quantity and we do not change the 
topology of the tube, this energy does not change. For 
the mean curvature, standard geometry consideration 10 
of the tube surface gives: 



H 



l + R'jz) 2 ~R{z)R"(z) 
2R(z)[l + R'{z) 2 ] 3 / 2 ' 



(8) 



Substituting Eq. JHJ| into the expression for bending en- 
ergy (the second term in Eq. (Q), and keeping terms up 
to second order in e, one obtains for the bending energy 
change: 



AE h 



e 2 ^-[3 + 2(kR ) i -(kRo) 



(9) 



III. ELECTROSTATIC ENERGY AND 
RELAXATION OF CHARGED SURFACTANTS 

Let us now proceed to calculate the electrostatic en- 
ergy change in the system under pearling deformation. 
We use the standard Debye-Huckel (DH) approximation 
to describe interactions between the charged surfactant 
molecules. In this approximation, the only role of free 
ions in solution is to screen the Coulomb potential of a 
charged surfactant. In other words, the electrostatic po- 
tential due to a charged surfactant molecule at a distance 
r from it is 



V UH (r) = eexp(-Kr)/Dr, 



(10) 



where e is the charge of one surfactant molecule (without 
loss of generality, we assume the charge of the surfactant 
is positive, D = 80 is the dielectric constant of water), 
and k is the inverse Debye-Hiickel screening radius. If 
the concentration of monovalent ions is water solution is 
cq, k is given by: 



^%-KCve 2 lDk B T 



(11) 
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The deformation of the tube also leads to the redistri- 
bution of charged molecules (or charge relaxation). The 
degree of charge relaxation depends on the specific sys- 
tem. In this section, for simplicity we assume the re- 
laxation of charged surfactant happens instantly and the 
distribution of surface charge is the equilibrium distribu- 
tion with respect to a given shape of the tube. We return 
to this assumption in detail in the next section. 

To find the equilibrium charge distribution, cr(z), 
which, in turn, enables us to calculate the change in the 
electrostatic energy, one needs to solve the DH equation 
for the electrostatic potential, V(r), in the system 



V 2 V(r) = n 2 V{v) 



(12) 



self-consistently with the boundary condition that the 
electric field at the tube surface is 2tt(t(z) and the sur- 
face charge is at a constant potential. In this paper, we 
use a simpler approach. Namely, we use a variational 
approach to calculate the electrostatic energy. Assuming 
the following ansatz for the charge distribution of the 
deformed tube: 



a(z) = c[l + xcos(kz)], 



(13) 



we optimize the electrostatic energy of the tube with re- 
spect to the variational parameter x. As we shall see 
later, x is proportional to e. This, coupled with the fact 
that for small deformation the response of the system is 
linear, the charge distribution obtained using variational 
approach is actually the true charge density of the system 
up to the second order in e. 

The conservation of the total charge of the tube: 

a 27ri? = — / / dS(z, <f>)a(z) , (14) 

27F J-TT/kJO 

immediately gives for the average charge density a: 
1 - (kR Q ) 2 



a ~ (T i 1 



X£ 



(15) 



The electrostatic energy of the tube with surface 
charge density, Eq. ljP3)l. is: 



Ee = 2 / rf ^( Zl '^ 1 ) rf ^( Z2 '^ 2 ) a ( Zl ) cr ( 2: 2) x 
V DH [d{ Zl ,z 2 ,fa,h)] 



(16) 



where the distance d(z%, z 2 , </>i> fa) between the two 
points (zi,fa) and (z2,fa) on the tube surface is 

d 2 (z 1 ,z 2 ,fa,fa) = [zi - z 2 ) 2 + R 2 {z 1 ) + 

R 2 (z 2 ) - 2R{ Zl )R(z 2 ) cos(fa - fa) (17) 



Substituting Eq. (|5I10I13I15|) and JUJ into Eq. lfHS)l. 
and expanding the integrand to second order in e (x and 
£ are of the same order of smallness), one obtains the 



following expression for the electrostatic energy per unit 
length of the tube after integration: 



2tt(j 2 R 2 

E e ~ — ^— ^(a + £ 2 a + £bx + cx 2 ), 

where the coefficients ao, a, b, and c are: 
a = 2ttIo(kRo)K (kRo), 



(18) 



nl (Vn 2 + k 2 R Q ) K [V^ + ^R, 



R 2 a d 2 3i? d 



4 dR 2 



4 dR 



«o 
2 



^kkR { 



°G 2 Uk 2 R 2 



1 

11-1 

2 ' 2 ' 2 



VAk 2 



^ kR °G 2 1 { k 2 R 2 



(k 2 + K 2 )R 2 





-1 1 -1 

2 ' 2' 2 



V<k~ 2 



K ) R C 21 



b = 2c + R ( 



2 

dRo \ 2 
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(fc 2 + k 2 )R 2 



1 

11-1 
2 ' 2 ' 2 





-1 1 -1 
2 ' 2 ' 2 



(19) 



Io and Kq are the modified Bessel functions of zeroth 

^' ~\ is the Meijer's G functioni^. 

Minimizing the electrostatic energy, Eq. I)18[l. with 
respect to x, one gets for x and the electrostatic energy 
change per unit length: 



x — ~£b/2c. 

D 



a Q = £ 



D 



4c 



(20) 



As expected, x is of the same order of smallness as e. 
This is consistent with the starting assumption we use in 
the expansion, Eq. (fT%|l . 

If the charged surfactant molecules do not relax to 
equilibrium surface distribution, their density remains 
constant under the deformation, x — 0. From Eq. i|18|l . 
the change in the electrostatic energy in this case is given 
by: 



= £ 



2 2k<j 2 R 2 



D 



(21) 



Correspondingly, the energy change due to the relaxation 
of charged surfactants comes from the two x-dependent 
terms in Eq. (|18[1: 



_ 2 2na 2 R 2 b 2 
D 4c 



(22) 



As one sees from Eq. (I19II . c is a positive coefficient. 
Thus the relaxation energy is negative as expected: elec- 
trostatic relaxations lower the cost of deforming a tube 
into pearls. 
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IV. DISCUSSION 

In this section, we comment on the relative importance 
of various energies in the system starting with the elec- 
trostatics energy and the contribution coming from the 
relaxation of charge molecules given by Eq. I|20[l and Eq. 
(|22[) . To gain a better physical insight into these equa- 
tions, it is instructive to consider the strong screening 
case, k ^S> k and kRo ^> 1 and expand the energies in 
powers of 1/kRq. For Ai?" orcl , the zeroth order term of 
the expansion is 



^^norcl (0) 



2D 



^ [(kR f - 1] e\ 



(23) 



It is easy to see that this energy behaves in the same way 
as the surface tension energy, Eq. Q. One, therefore, 
identifies the electrostatic contribution to the surface ten- 
sion of the surfactant tube: 



l /D. 



(24) 



It is not surprising to see that, in the absolute value, this 
"electrostatic" surface tension is simply the electrostatic 
energy per unit area of a flat bilayer at the same charge 
density. The negative sign in this expression reflects the 
fact that electrostatic repulsions between charged surfac- 
tants prefer to increase the area of the surfactant bilayer. 
The total surface tension of the layer, of course, remains 
positive because of the non-electrostatic interaction be- 
tween surfactant molecules counter-balance this negative 
electrostatic contribution. 

The next non-zero term of the expansion of AE^° lcl is 
of second order in 1/kRq: 



a pnorel (1) r 

e 

{[2~2(kR ) 2 



2 2 
7T (if. hZ 



_1 i?0 



f 



16£> (kR ) 2 
- [3 + 2(fci? ) 4 - 



(fci?o) 2 ]}(25) 



The first square bracket term simply adds a small cor- 
rection, a e / 4(nR ) 2 , to the electrostatic surface tension 
energy. Comparing the second square bracket term with 
the bending energy change, Eq. one immediately 

identifies this term as the "electrostatic" bending energy 
change with the corresponding "electrostatic" bending 
rigidity given by 



7c 



J /4D. 



(26) 



It is positive and, within a constant numerical factor, 
agrees with the well-known expression 15 for 7 C calculated 
using other methods. Thus expanding the electrostatic 
energy with respect to the tube curvature (1/kRq) in 
this near flat limit, one recovers all standard formulae 
for the "electrostatic" contributions to the elastic param- 
eters (surface tension and bending rigidity) of the bilayer 
surface. 

Also in this limit, the relaxation energy becomes, to 
the lowest order in 1/kRq: 



ae: 



,Tra e Ro {kR ) 4 
16 (kR ) 4 



(27) 



Thus, the relaxation of charged molecules belong to 
the fourth order or higher in the expansion with respect 
to the tube curvature. Since the electrostatic surface 
tension and the electrostatic bending energy are, corre- 
spondingly, the zeroth and second order terms in this 
expansion (the first order term in the expansion vanishes 
because of the symmetry of the reference flat surface), 
charge relaxation energy gain is parametrically small 
compared to the electrostatic surface tension and bending 
energy in this limit and can be ignored. In other words, 
in this strong screening limit, the charge density of the 
surfactant bilayer can be considered uniform during the 
deformation of the tube. 

In the opposite limit of very weak screening, kRq <C 1, 
generally speaking, the electrostatic interaction is long 
range and is so large that linear analysis becomes in- 
valid in a very short time after the instability develops 
and nonlinear terms must be included in describing the 
development of instability. This is, however, a very com- 
plicated task. This is why in literature one usually as- 
sumes the final (pearl-on-a-string) structure of instability 
is given and variationally minimized its total energy to 
find its parameters (size, period). Nevertheless, for the 
discussion of the role of charge relaxation energy, one 
can still use the result a linear analysis given by exact 
expression, Eq. I)2()[l. which is valid at a very early time 
of instability. 

Expanding the energies in powers of kRo, for k > k, 
we get to the lowest order in kRq, 



AE t 
AE, 



_ £2 2TT 2 a 2 R 2 Io(kR ) [K (kRo) - kRoK^kRo)]' 



D 
D 



(kR ) 2 



h(kSo)Ki(kRo) - ^ 



Ko(kRo) 
I {kRo)K (kR ) 
Io{kR ) 



K (kR 
Ko(kR Q ) — kR Q K 1 (kR Q ) J j . 



-kRoK^kRo) x 



(28) 



Since for A; > n, all length scales are smaller than the 
screening radius and electrostatic interactions are not 
screened. Correspondingly, the energies are independent 
of k as shown by the above equation Q28JI. 

For k <C k, we get to the second lowest order in kRq 
and k/n: 



AET = 



£2 2^| ln ^o +(2 + 7£)+ 



D 

IP 

2 k 2 k 2 \tl(k,Ro/2) J 



IE 



AE, = e 



2 2n 2 a 2 Rl r 



D 



1 - 



IE 



k 2 \u(kR /2) 



(29) 



Here = 0.5772 is the Euler's constant. For the later 
case, we see that the relaxation energy is larger than the 
total electrostatic energy by a large logarithmic term, 
ln(fti?o/2). This is because, the first expansion term 
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which logarithmically diverges with kRq — > in the re- 
laxation energy is exactly equal in magnitude and op- 
posite in sign to the the first expansion term for the 
non-relaxation energy. As a result, the total electrostatic 
energy, contains only the second and higher order expan- 
sion terms. Thus, it is parametrically smaller than either 
of these components. Obviously, within linear analysis, 
the inclusion of the relaxation energy is important in this 
limit to get the correct behaviour of the electrostatic en- 
ergy. 

For the experimental situation of Ref. where the 
screening radius, the tube radius and the pearl size are 
comparable to each other, the situation is marginal. 
Therefore, one might ask whether or not the charge relax- 
ation still plays a significant role. To answer this ques- 
tion, we numerically evaluate the exact (within linear 
analysis) expressions, Eqs. lj?T]l and for the A_E" oro1 
and A_Eg Cl , respectively. In Fig. [31 we plot the ratio be- 
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FIG. 3: The absolute value of the ratio between the electro- 
static energy gain due to the relaxation of charge molecules 
to the total electrostatic energy when a tube deforms into a 
pearling structure as a function of the wavevector of defor- 
mation, kRo for kRo = 1. Four different values of kRq — 0.1, 
0.5, 1.0 and 2.0, are used. Lighter curve corresponds to higher 
Ki?o- The divergencies observed near kRo ~ 1 is due to the 
vanishing of AE e owing to Raleigh instability. 

tween the energy gained due to the relaxation of charged 
molecules and the total electrostatic energy change for 
different values of kRq. The divergence of this ratio 
observed at about kRo ~ 1 is due to fact that the to- 
tal electrostatic energy change goes through zero at this 
wavevector. As one can see from this figure, for strong 
screening kRq > 1, the relaxation energy contributes a 
small part of the total electrostatic energy while for weak 
screening kRq < 1, it contributes significantly to the total 
energy. However, for the marginal case kRq = 1 which 
is more relevant to the experiments, the numerical re- 
sult shows that the relaxation of charged molecules con- 
tributes about 20%, which is a small fraction. For smaller 
screening radius, the exponent of 4 in Eq. (|27|l quickly 



reduces the relaxation energy to an irrelevant contribu- 
tion in the total energy. Thus, one concludes that in the 
experimental situation of Ref. H the electrostatic relax- 
ation is a small effect and the pearling of the tubes is 
dominantly due to the Raleigh instability. 

Before concluding, let us come back to the assumption 
made in the previous section that the surface charge is 
always in equilibrium with a given shape of the tube. 
To show that this is a reasonable assumption, let us es- 
timate the charge relaxation of the surfactant bilayer. 
Indeed, one can view the tube as an RC circuit. One 
typical wavelength ~ Rq 7 the conductance of the cir- 
cuit is the same as the conductivity, TlT 1 ~ aoe/j, with 
fi ~ 10 9 m/sN the mobility of the surfactant. Because all 
charges are screened at the distance the capacitance 
of this circuit is of the order C ~ R 2 k. Thus the relax- 
ation time of this circuit (or of our surfactant tube) is of 
the order (TZC)^ 1 ~ <7oe/x/i? 2 K which is about 0.1ms us- 
ing experimental parameters. This is much smaller than 
the growth rate of the pearling instability (in seconds). 
Thus the surfactant charges, to a good approximation, 
can be considered always in equilibrium. 



V. CONCLUSION 

We have considered the relaxation of charged surfac- 
tant in the pearling instability for surfactant tubes ob- 
served in the experiments of Ref. 0. This effect is known 
to be the main driving force for pearling instability of 
polyelectrolytcs in poor solvent for which screening is 
very weak. Using linear analysis, we showed that for the 
marginal situation of Ref. |5|, the effect is small and the 
pearling instability is mainly due the Raleigh instability 
caused by the surface tension of the bilayer. 

Charge relaxation becomes important only for weak 
screening such that kRq <C 1. In this limit, in our linear 
analysis, the relaxation energy is larger than the total 
electrostatic energy by a logarithmic factor, Ih(kRq/2). 
For very weak screening condition, its inclusion is a must 
if one wants to obtain the correct behaviour of electro- 
static energy. Thus, although charge relaxation plays 
small role in the biological phenomena of Ref. 0J, it might 
well play a role in the domain of flowing microemulsions, 
where the tube diameters are much smaller. 
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